Perturbation theory of the mass enhancement for a polaron coupled to acoustic 
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We use both a perturbative Green's function analysis and standard perturbative quantum me- 
chanics to calculate the decrease in energy and the effective mass for an electron interacting with 
acoustic phonons. The interaction is between the difference in lattice displacements for neighbour- 
ing ions, and the hopping amplitude for an electron between those two sites. The calculations are 
performed in one, two, and three dimensions, and comparisons are made with results from other 
electron-phonon models. We also compute the spectral function and quasiparticle residue, as a 
function of characteristic phonon frequency. There are strong indications that this model is al- 
ways polaronic on one dimension, where an unusual relation between the effective mass and the 
quasiparticle residue is also found. 



I. INTRODUCTION 



When electrons interact strongly with phonons, the 
electrons acquire a polaronic character, i.e. they move 
around the lattice much more sluggishly than non- 
interacting electrons would, because a polarization cloud 
must accompany them as they move. A measure of the 
strength of the coupling between the electron and the 
phonons is the degree to which the ground state energy 
is lowered. For example, previous studies for the Hol- 
stein modeli have indicated that the decrease in energy is 
proportional to the bare coupling strength (A) in strong 
coupling, 2 independent of the value of the phonon fre- 
quency. On the other hand, in weak coupling, while the 
proportionality to A remains, there is some dependence 
on phonon frequency, and in fact, the decrease in energy 
is greater for higher phonon frequency!^ 

A much more indicative measure of the polaronic char- 
acter of an electron is the effective mass. In the Holstein 
model a glimpse of polaronic tendencies, even within per- 
turbation theory, can be attained by examining the ef- 
fective mass, particularly in one dimension. Usually an 
increasing effective mass is accompanied by a decrease 
in quasiparticle residue, although this is not always the 
case, as described below. 

The Holstein model describes electrons interacting 
with optical phonons; the coupling is via the electron 
charge density, and, in this sense, the Holstein model 
serves as a paradigm for electron-phonon interactions just 
like the celebrated Hubbard model^ is the simplest de- 
scription of electron-electron interactions. Many of the 
basic features of this model are now fairly well under- 
stood — see Ref. @and@ along with more recent work 
in Ref. [H andQ- However, just as important is the 
electron interaction with acoustic phonons; typically the 
ionic motions couple to the electron motion, as opposed 
to its charge density. A very simple model to describe 
this kind of electron-phonon interaction within a tight- 
binding framework is given by 



(i) 



where angular brackets denote nearest neighbours only, 
and 
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)<5i,i±a x - a(u yi - u yj )Sij±a y ■ (2) 



This Hamiltonian has been written specifically for two di- 
mensions, but the generalization to three dimensions (or 
back to one dimension) is evident from Eqs. (JlJ and (|2]l. 



The operators and parameters are as follows: 



creates (annihilates) an electron at site i with spin a. 
The x-components for the ion momentum and displace- 
ment are given by p X i, and displacement u x i, respectively 
(similarly for the y-components) , and the ions have mass 
M and spring constant K connecting nearest neighbours 
only. The electron-ion coupling is linearized in the com- 
ponents of the displacement, and we choose to include 
only longitudinal coupling. 

This Hamiltonian is commonly known as the Su- 
Schrieffer-Heeger (SSH) model^^ because it was used 
for seminal work describing excitations in polyacetylene 
by these authors. However, it was also introduced and 
studied a decade earlier by Barisic, Labbe, and Fricdcl 10 
to describe superconductivity in transition metals, so 
we will refer to it as the BLF-SSH model. Much of 
the work done on this model is in the adiabatic ap- 
proximation, i.e. the phonons are treated classically. 8,9 
This was followed by an examination of quantum fluctua- 
tions through quantum Monte Carlo and renormalization 
group studies^ and these authors focused on half-filling. 
They found that the lattice ordering (in one dimension) 
was reduced by quantum fluctuations. 

Very little work has been done, however, in the quan- 
tum regime for a single electron. Capone and coworkers 
studied a model similar to this one, except that they uti- 
lized optical phonons instead of acoustic ones j 12 ' 13 This 
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leads to some significant differences, about which we will 
comment below. In the past decade Zoli has studied 
the BLF-SSH polaron using perturbation theory, and 
found, for example, a perturbative regime in one dimen- 
sion where polaron effects are absent^ This result hap- 
pened to agree with the conclusions of Capone et alJ^ 
in the perturbative regime of the CSG models In this 
paper we focus on 2nd order perturbation theory, and 
find results in disagreement with Ref. (TiJ . These results 
also disagree qualitatively with the results from the CSG 
model. That is, in one dimension, for example, pertur- 
bation theory breaks down as the characteristic phonon 
frequency decreases. In two dimensions there is a modest 
mass enhancement for all characteristic phonon frequen- 
cies, while in three dimensions the mass enhancement ap- 
proaches unity in the adiabatic limit. We also note that 
the quasiparticle residue does not necessarily follow the 
trend of the inverse effective mass, as the characteristic 
phonon frequency varies. 

This paper is organized as follows. In the following sec- 
tion we outline the calculation, both using perturbation 
theory, and using Green function techniques. For some 
of our work (especially in one dimension), the calcula- 
tion can be done analytically, and we derive these results 
where applicable. In Section III we show some numeri- 
cal results and compare our results with previous work 
and other electron-phonon models. We close in the final 
section with a summary. The main conclusion is that, 
as far as one can tell from weak coupling perturbation 
theory, the BLF-SSH model has a stronger tendency to 
form a polaronic state than is the case with the Holstein 
model. In one dimension this is most evident in the ef- 
fective mass, and not at all evident in the quasiparticle 
residue. 



II. PERTURBATION THEORY 

A. Hamiltonian 

The Hamiltonian Eq. ([IJ, Fourier-transformed to 
wavevector space, and utilizing phonon creation and an- 
nihilation operators, is written (again in 2D), 

ka 

a xq a xq + a yq a yq] 
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Here, 

e k = e(k x , k y ) = -2i[cos (k x ) + cos (fc y )] (4) 



is the dispersion relation for non-interacting electrons 
with nearest neighbour hopping, and 
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is the phonon dispersion for acoustic phonons with near- 
est neighbour spring constants K, and ujq = \/4K/M 
is the characteristic phonon frequency. The phonon cre- 
ation and annihilation operators are given by a xq and a xq , 
respectively, and similarly for those in the y-direction. 
The coupling "constants" are given by 



g x (k, k') = ia. 



MNui(k - ft') 



sin (k' x ) — sin (k x ) 



(6) 



with a similar expression for the y direction, and M is 
the mass of the ion and N is the number of lattice sites. 



B. Green's function analysis 

Carrying out a Green's function analysis using the free 
electron and phonon parts of the Hamitonian as the un- 
perturbed part, gives, for the self energy of a single elec- 
tron to lowest (2nd) order in the coupling a, 

E(ft, uj + iS) = 

-X>*(*> k')\ 2 + \g y (k, k')\ 2 ]G (k',u J + iS ~ tj(k - k')), 



(7) 



where Go(k,uj + iS) = \u + iS — e^l is the non- 
interacting electron retarded propagator. 

One way to determine the effect of interactions on the 
electron dispersion is to compute the renormalized en- 
ergy for the ground state (here, k x = k y = 0), and the 
effective mass. The effective mass has long been used 
as the primary indicator for polaronic behaviour^*, and 
though within 2nd order perturbation we can only get 
an indication of this crossover, we use it here nonethe- 
less. The renormalized energy is given by the solution 
for the pole location in the interacting electron Green's 
function, G(k, cu + iS) — \ui + id — e k — S(ft, w + iS)] , 



E k = e k + ReX(k,E k ). 



(8) 



To determine the effective mass, defined by the expecta- 
tion that Ek = h 2 k 2 /(2m*), we take two derivatives^ of 
Eq. ©, and, using the fact that (dE k / dk)\k=o — 0, we 
obtain 



m 1 I i d 2 s(fc,") I 

dZ(k,uj) 1 d 2 Z(k,uj) 

dej |w=Efc " 2t dk 2 ^ =Ek 



Here we have used the fact that the band mass given by 
the electron dispersion in Eq. Q is to = l/(2t). Note 
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that it is common (and advisable) to replace the substitu- 
tions for ui required in Eq. Q with e k , rather than with 
Ek- This is due to the fact that the former substitution 
keeps the evaluation for every term at 0(a 2 ), whereas the 
latter substitution includes some (inconsistently) higher 
order contributions. The former substitution is known as 
Rayleigh-Schrodinger perturbation theory while the lat- 
ter is known as Brillouin-Wigner perturbation theory*^ 
This means that we will use the following equation, 



that this corresponds to Rayleigh-Schrodinger perturba- 
tion theory with the self energy, evaluated at u) = e k cor- 
responding to the 2nd order energy correction. Eq. (fT"3|) 
can be evaluated numerically, and then two derivatives 
with respect to k are required. However, the same nu- 
merical problems mentioned above will arise; fortunately, 
at least in one dimension, Eq. (| 13[) can be evaluated an- 
alytically, whereas we were unable to do the same with 
Eq. ©. 
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1 d 2 Z(k,uj) 
2b dk 2 1 



(10) 



to define the effective mass. 

In contrast the quasiparticle residue is defined as the 
weight that remains in the (5-function-like portion of the 
spectral weight. The spectral weight is defined as 



A(k, uj) = --ImG(fc, uj + iS) 

7T 

-Im- 



7r uj + i8 — efe — S(fc, uj + iS) 



(11) 



For a given momentum, as the energy of the pole given 
by Eq. ([5]) is approached, the imaginary part of the self 
energy tends towards zero; this produces a ^-function 
contribution in Eq. (11 1 [) , at the pole energy, but with 
weight Zk defined by 



Zk 



1 



1 



d£(fc,u) | 



(12) 



The relationship amongst these various quantities — ef- 
fective mass in Eq. ©, effective mass in Eq. ((TOl) . and 
quasiparticle residue in Eq. (fT2|) — is discussed further 
in the Appendix. 



C. Standard perturbation theory 

Eq. ([9]) requires a numerical evaluation of Eq. (J7]), 
and then the required derivatives can be (numerically) 
determined. Because the positions of the singularities 
in Eq. ([7]) are difficult to determine in advance, it is 
customary to introduce a small (numerical) imaginary 
part corresponding to the infinitesimal <5, and then the 
numerical integration is more stable. This trick remains 
problematic, as we discuss further below. Alternatively, 
we can simply perform a 2nd order perturbation theory 
expansion, as outlined in every undergraduate quantum 
mechanics textbook. The result is 
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(2) = 

k M N ^ 

k' 



uj(k-k') [e k - e k > - u>(k - k 1 )] 



(13) 

where we remember that the first order (in a) contribu- 
tion is of course zero, and the superscript (2) indicates the 
2nd order contribution. Comparison with Eq. Q shows 



III. RESULTS AND DISCUSSION 
A. Analytical results in ID 

The result of an analytical evaluation^ of Eq. (|13p is, 
in one dimension, 

£ (2) (fc) = -— AblfWo{-2cosA; + 7rDci + C k (ui Q )\, 



(14) 

where (Do = u>q/ (At) , and a dimensionless coupling pa- 
rameter Ablf is defined, in analogy to the dimensionless 
coupling parameter defined in the Holstein model, as 



A 



BLF 



a 2 1 

Mui 2 W' 



(15) 



where here the bandwidth W — At for one dimension. 
Note that this coupling parameter has nothing to do 
physically with the coupling parameter defined in the 
Holstein model, so we will treat them as completely 
independent.— The function Ck(uio) must be evaluated 
separately in the two regimes: 

C fc (w ) = 2^Juj 2 - lfh(k)+h(-k)-2h{-K/2)j, w > 1, 

(16) 

where 



^)=tan-^"° t / a ^ + 1 

V v^o-i 



(17) 



and 



C k {w ) = ^l-Co 2 (s{k) + s{-k) - 2s(tt/2)^ , Co < 1, 

(18) 

where 



s(k) = lo| 



/ cj tanf + 1 + yjl-wl 
'yw tan| -I- 1 - yr 



r,2 



(19) 



Eq. ([141) is readily evaluated at k — to determine the 
ground state energy. Evaluating the second derivative 
with respect to wave vector k is equally straightforward, 
and determination at k = yields the rather simple result 
for the effective mass, 

32 Ablf 



m 



= 1 



TO 7T U>0 

valid for all values of CjqJ^ 



(20) 
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Comparison with other models 



An analytical result is readily available for the Holstein 
model; there, the ground state energy (in ID) was given 
by2 



E H = -2t[ 1 + A H 



uj e 



Cje + 1 



(21) 



where uje = U)e/(4£) is the Einstein phonon frequency 
normalized to the bandwidth, and, as explained earlier, 
the dimensionless coupling constant Ah cannot be com- 
pared directly to the corresponding quantity for the BLF- 
SSH model. The effective mass is given by 



in" 
m 



1 + 2Cj e 




(22) 



In both cases, as the characteristic phonon frequency ap- 
proaches zero (adiabatic limit) the ground state energy 
approaches the non-interacting value; however, the effec- 
tive mass diverges in this same limit. So, while the first 
statement would appear to justify perturbation theory in 
this limit, the second statement clearly indicates a break- 
down in the adiabatic limit. It is known in both cases 
that the adiabatic approximation leads to a polaron-like 
solution for all coupling constants ; 20 ^ 1 and clearly these 
two observations are consistent with one another. In fact, 
the divergence is stronger in the BLF-SSH model, and 
goes beyond the inverse square-root behaviour observed 
for the Holstein model and attributed to the diverging 
electron density of states in one dimension^ this indi- 
cates that the BLF-SSH model, at least in the adiabatic 
limit in one dimension, has a stronger tendency for po- 
laron formation than the Holstein model. 

Interestingly, in the model studied by Capone et alM, 
where optical phonons were used, the opposite behaviour 
was obtained; they found that the effective mass ratio ap- 
proached unity as the characteristic phonon energy ap- 
proached zero.— In the opposite limit Capone et alj^ 
found an effective mass ratio that did not approach unity 
as the characteristic phonon frequency increased (anti- 
adiabatic limit). In the BLF-SSH model, however, this 
ratio does approach unity as the phonon frequency in- 
creases beyond the electron bandwidth, in one dimen- 
sion, in agreement with the Holstein result in all dimen- 
sions. As we will see below, however, in the BLF-SSH 
model in two and three dimensions the effective mass ra- 
tio remains above unity in the anti-adiabatic limit. This 
is not surprising, since here the interaction modulates 
the hopping, and we expect a non-zero correction in this 
limits In the adiabatic limit, the BLF-SSH mass ratio 
approaches a constant value in two dimensions, and falls 
to unity in three dimensions, both in agreement with the 
behaviour in the Holstein model. 

Our results disagree with those of ZoliM for reasons 
that are not entirely clear. We have utilized both the 
straightforward perturbation theory method (analyti- 
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FIG. 1. Electron self energy for the ground state (k = 0), 
normalized to A (or Ah) vs. characteristic phonon frequency 
(jjo (this is uje for the Holstein model), for both the BLF-SSH 
and Holstein models, in one, two, and three dimensions, as in- 
dicated. Alternatively, the ordinate is simply the second order 
(in g) correction to the ground state energy within Rayleigh- 
Schrodinger perturbation theory. In all cases the magnitude 
of the correction increases with increasing uiq. At low uiq the 
magnitudes of the the results are ordered 3D, 2D, ID (low- 
est to highest) whereas at high frequency the ordering is just 
the opposite. All six cases have non-zero limiting values as 
ui{_) — > do, given in Table 1. 



cally and numerically), and the Green's function formal- 
ism (numerically). In the latter case we required a nu- 
merically small imaginary part for the frequency signif- 
icantly smaller than the value quoted in Ref. (|14l ) (we 
used S = 10~ 9 whereas he used 5 = 10~ 4 . However, 
as is clear from our analytical result, Eq. (|2"0)) . our ef- 
fective mass diverges at low phonon frequency, and de- 
creases monotonically to unity as the phonon frequency 
increases. The result in Ref. ((Til ) peaks sharply near 
u)q f=3 1, and, as noted above, decreases to unity at low 
phonon frequency. 



C. Numerical results 

In Fig. Q] we plot the reduction in the ground state 
energy due to the second order correction (for the BLF- 
SSH model, this is given by Eq. (fT3|) L normalized to A 
(or Xh)- This is also written as E(fc = 0,oj = £fc)/A, 
where the self energy is given by the expression in Eq. 
([7]). Also plotted for comparison are the corresponding 
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quantities for the Holstein model. Note that both models 
share a few features in common: (i) they both go to zero 
as the characteristic phonon energy decreases to zero, re- 
gardless of the dimensionality, (ii) they all approach a 
non-zero negative (and finite) value as the characteris- 
tic phonon frequency grows, and (iii) they cross one an- 
other in strength as a function of dimensionality as ujq 
increases, i.e. at low phonon frequencies the self energy 
has the highest magnitude for one dimension, whereas for 
high phonon frequency the highest magnitude is achieved 
in both models for three dimensional systems. Also note 
that the BLF-SSH results are well separated from Hol- 
stein results. In particular, there appears to be more 
'bang for the buck' with the BLF-SSH model, i.e. for a 
given value of A and the same characteristic phonon fre- 
quency, the energy reduction is almost an order of mag- 
nitude higher for the BLF-SSH model as compared with 
the Holstein model. Again, we remind the reader that the 
value of A in the Holstein model has nothing to do with 
the value of A in the BLF-SSH model, so this comparison 
is unwarranted. 

For this reason we will use the value for the self energy, 
in weak coupling, as the phonon frequency increases to 
infinity, as the energy scale that provides a measure of 
the energy lowering expected for a given model and a 
given dimensionality. These numbers, mostly determined 
analytically, are provided in Table I. 

TABLE I. lim^oo E(fc = 0,w = e k )/(Xt) 
Dim. BLF-SSH Holstein 



ID 


-16 


-2 


2D 


-23.3 


-4 


3D 


-30.2 


-6 



In Fig. [5] we plot the effective mass ratio (minus unity), 
normalized to the self energy evaluated for infinite char- 
acteristic phonon frequency. This normalization is im- 
portant to divide out enhancements that are solely due 
to definitions. Moreover, in this way, we are determining 
the mass enhancement for a given 'coupling strength', 
where this strength is now a measure of the energy low- 
ering caused by a certain amount of coupling to phonons, 
regardless of the origin of that coupling. This plot now 
makes clear that the BLF-SSH model, within weak cou- 
pling perturbation theory, has more 'polaronic' tendency 
than the Holstein model. Note in particular that the di- 
vergence (in ID) at low characteristic phonon frequency 
is much stronger for the BLF-SSH model, as Eq. (pTO)) 
already indicated. Thus, as discussed above, we antici- 
pate that in the adiabatic approximation, in ID, the sys- 
tem will always be polaronic, regardless of the coupling 
strength, in agreement with the result of the Holstein 
model^ and in disagreement with the result from the 
hybrid model defined in Ref. HH. 

Otherwise, the behaviour of the effective mass in the 
two models is very similar, as a function of characteris- 




FIG. 2. The electron effective mass, normalized to the 2nd 
order correction to the energy for the anti-adiabatic limit, vs. 
characteristic phonon frequency, wo, for both the BLF-SSH 
and Holstein models, in one, two, and three dimensions, as 
indicated. In ID the effective mass diverges for both models, 
though the divergence is stronger for the BLF-SSH model, as 
indicated by Eq. (|20[1 . In 2D the effective mass approaches a 
constant as ujq — > for both models, while in 3D the effective 
mass ratio approaches unity in the same limit. At the opposite 
extreme, both ID results give m*/m — > 1 as ujq — > oo, while 
in both 2D and 3D the effective mass remains above unity 
in this limit. Note that in all three dimensions, for a given 
reduction in energy as given by the 2nd order correction to 
the energy, the BLF-SSH model results in significantly higher 
effective masses. 

tic phonon frequency, for the various dimensions shown. 
The effective mass can be made arbitrarily close to unity, 
for any non-zero phonon frequency, for sufficiently weak 
coupling. Preliminary numerical calculations indicate a 
free electron-like to polaron crossover,— similar to what 
was found for the Holstein model. 



D. Spectral function 

It is interesting to examine the spectral function, de- 
fined by Eq. (TlTT) (see also the discussion in the Ap- 
pendix). For simplicity we show the result in one dimen- 
sion, in Fig. 02 for the ground state (k = 0) as a function 
of frequency. 

The results for two or three dimensions do not differ 
in any significant way from these results. The results 
for three different characteristic phonon frequencies are 
shown. In each case a quasiparticle ^-function is present 
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FIG. 3. Spectral function for the BLF-SSH model, for A = 0.2 
for three different characteristic phonon frequencies, as a func- 
tion of frequency. All three spectra are similar as one would 
find for the Holstein model, and consist of quasiparticle peak 
with weight z = 0.766,0.727,0.724, for u /t = 0.1,0.5,2.0, 
respectively, followed by an incoherent piece. 




FIG. 4. Quasiparticle residue, zo vs. coo/t for both the BLF- 
SSH and Holstein models, in all three dimensions. Note that 
while the result for the Holstein model tends to be inversely 
proportional to the effective mass, this is not the case for the 
BLF-SSH model at low phonon frequency, and in ID and 2D. 
In one dimension in particular, the effective mass diverges, 
while 2q also turns upward. 
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FIG. 5. Comparison of the quasiparticle residue (upper 
panel) with the electron effective mass (lower panel) as a func- 
tion of LJo/t, for the BLF-SSH model in one dimension. The 
behaviour noted in Fig. [4] is clear here. Moreover, note the 
scales; while the effective mass ratio is very large (~ 4) for 
A = 0.01 and small values of ujo/t, the quasiparticle residue 
remains within 15% of unity. 



(here artificially broadened so as to be visible), followed 
by an incoherent piece; the incoherent part has energies 
ranging approximately from — 2t < u> < +2t + uq. The 
quasiparticle residue, zq must be determined numerically, 
and is given in the figure caption for each of the cases con- 
sidered (see also Fig. 0}. We have verified that the re- 
maining weight (the spectral functions each have weight 
unity) is present in the incoherent part. The result shown 
is not too different from what is found in the Holstein 
model; the singularities from the ID electron density of 
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states are now smeared out in the incoherent piece, as 
a result of the coupling and phonon energy having some 
frequency dependence. We show in Fig. 21 as a function 
of wo, the quasiparticle residue for both the Holstein and 
BLF-SSH models. The Holstein results tend to follow the 
inverse of the result for the inverse effective mass; this is 
as expected. This is not the case with the BLF-SSH, but 
for more subtle reasons than the fact that the self energy 
is now momentum dependent. The more important ef- 
fect, which shows up in both ID and 2D results, is that 
the quasiparticle weight requires an evaluation of the fre- 
quency derivative of the self energy at the energy of the 
pole, whereas the effective mass in Raylcigh-Schrodinger 
perturbation theory requires the same derivative at the 
non-interacting ground state energy. Most noteworthy 
is that the quasiparticle residue shows a clear upturn at 
low characteristic phonon frequencies, while the inverse 
effective mass clearly approaches zero (see Fig. [5]) as this 
characteristic frequency is taken to zero. 

To see this more clearly we show in Fig. [5] a compari- 
son of the residue (upper panel) vs. effective mass (lower 
panel), as a function of oj , for two (weak) strengths of 
electron phonon coupling. At high phonon frequency, as 
the former decreases, the latter increases with decreas- 
ing phonon frequency, but at low phonon frequency, the 
two properties no longer behave in inverse fashion with 
respect to one another. 



IV. SUMMARY 



The BLF-SSH model appears to have very strong po- 
laronic tendencies, stronger than those of, say, the Hol- 
stein model, especially in one dimension. This conclusion 
is based on the 2nd order perturbative calculation per- 
formed in this paper, but also has corroborative evidence 
from calculations in the strong coupling regime. In one 
dimension we have been able to obtain an analytical so- 
lution for the ground state energy and the effective mass. 
The conclusion concerning polaronic behaviour is an im- 
portant one, as much of what we know about polarons 
arises from Holstein-like models^ In particular, for a 
coupling strength that leads to a fixed amount of energy 
lowering (in 2nd order), the effective mass can become 
an order of magnitude larger than the bare mass, a clear 
indicator that perturbation theory breaks down. This 
occurs in the BLF-SSH model at much weaker coupling 
than in the Holstein model. We have also noted that 
the relationship between effective mass and quasiparticle 
residue breaks down in one and two dimensions for the 
BLF-SSH model, not because of the momentum depen- 
dence in the self energy, but because the two properties 
involve evaluation of the frequency derivative of the self 
energy at different energies. Future work will address the 
strong coupling regime. 
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Appendix A: Perturbation Theory 

It is sometimes stated that for a momentum- 
independent self energy, the quasiparticle residue is equal 
to the inverse of the effective mass. This follows simply 
by comparing Eqs. © and (TT2"j) . On the other hand, we 
have argued that Eq. (fTU)) is more appropriate for the ef- 
fective mass, in which case this statement appears not to 
be true. A resolution of this difficulty is straightforward 
for the Holstein model, which we outline below, but, in- 
terestingly, not possible for the BLF-SSH model, at least 
in one dimension. The essential difference appears to 
be that in the Holstein model the (phonon) excitations 
are gapped, whereas they are not in the BLF-SSH model 
because of the low-lying acoustic modes at small momen- 
tum transfer. In this appendix we focus attention on one 
dimension, where some subtleties arise. 

For the Holstein model the computation of the self en- 
ergy in weak coupling is straightforward^ We obtain 



we 



(2ty 



(Al) 



The location of the quasiparticle pole at zero momentum 
(ground state) is then given by 



uj + 2t 



v / (w _ WB )2_ (2i)2 ' 



(A2) 



which can readily be determined numerically. Denoting 
the solution by writing uj = — 2t — Eb (so Eb is the 'bind- 
ing' energy below the bottom of the band), we can then 
use this in the spectral function, Eq. (fTTj) . to determine 
the residue z in the quasiparticle peak at uj = — It — Eb- 

A(k = 0, uj) = zq5(uj + 2t + Eb) + incoherent part. (A3) 

Straightforward calculation gives 

/ 2\ H uj E [l + 2Q E + 2Eb\ \ 

20 = 1/1 + - ; r^W . A4 

V [(I + 2^ + 2£ b ) 2 -l] 3/2 / 

which is not in agreement with the inverse of Eq. (|22[) , 
except when \h is truly very small. Here Eb = Eb/(4t). 

In particular, for arbitrarily small Xh, 
dY,{uj) / duj\ u =_2 t , which is used in Eq. (|2"2"|) . di- 
verges as uje — > 0, leading to a divergent effective mass 
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(and therefore associated residue of zero). On the other 
hand, from Eq. (|A2[) one readily sees 



lim E b = t(Xuj E /t) 



2/3 



from which Eq. (IA4|) yields the result 

lim zq — 2/3, 
(jj e — >o 



(A5) 



(A6) 



surprisingly a universal number. The actual weight in the 
quasiparticle peak of the spectral function given by Eq. 
pip for any given (even very small) value of \h actually 
tracks Eq. (|A4|) . and not the inverse of Eq. (|2"2"j) . 

Interestingly, for the Holstein model, one can take 
a different tact towards calculating the spectral func- 
tion: using perturbation theory to compute the per- 
turbed wave function, which is then inserted into the 
calculation for the matrix elements required in the defi- 
nition of the spectral function,— one obtains 



A pcrt {k = 0,Lj ) = zg ect 5(u + 2t + 



Xh^e 



A /(1 + 2£ B )2)-1 
1 2tuj E X H 6(2t- \w-u E \) 



tt(uj + 2t) 2 ^(2tf-{uj-u E f 



(A7) 



Note that there is no difficulty in integrating over this 
function, as the divergence in the denominator (l/(w + 
2t) 2 ) is not within (or bordering) the range of frequency 
given by the Heaviside function restriction in the numer- 
ator. This is due to the finite phonon frequency, loe- 



From this expression fulfillment of the sum rule deter- 
mines that 



pert 



= 1/1 + 



2\ H u) E [l + 2Q E ] 
[(1 + 2^)2-1] 3/2 



(A8) 



which is in agreement with the inverse of Eq. (J22J) . The 
message is that, as long as we use the expression given 
by Eq. ([TTjl for the spectral function, the area under the 
quasiparticle peak will correspond to Eq. (|A4[> . which 
is not the inverse of the effective mass, even if the self 
energy is independent of momentum. 

In the BLF-SSH model, the self energy is evaluated 
numerically through Eq. ((7]). An attempt to follow 
the procedure just outlined, which leads to Eqs. (IA7|) 
and (|A8|) for this model fails; this is because the mini- 
mum phonon frequency is zero, so the restriction corre- 
sponding to the Heaviside function in Eq. (IA7|) yields 
— 2t < lj < 2t + ujq] this in turn makes the divergence 
at lu — —2t non-integrable. One can only (in ID) define 
the spectral function through Eq. in which case 

the inverse of the effective mass differs from the quasi- 
particle pole for two reasons: the usual reason that the 
explicit momentum dependence now plays a role (see Eq. 
(fTUlO . and, in addition, the derivative of the self energy 
with respect to frequency is evaluated at u> = —2t for the 
effective mass, whereas it is evaluated at the frequency 
corresponding to the pole for the quasiparticle residue. 
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